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1. Introduction 

It is a paradigm of the AdS / CFT correspondence that the dynamics of a field theory living 
in an (asymptotically) anti-de Sitter bulk space-time encodes the correlation functions of 
its dual (deformed) conformal field theory jj], gj. The correspondence holds par excellence 
for J\f = 8 gauged supergravity in (d + 1) = 5 dimensions on the bulk side, which is 
obtained from D = 10 type IIB supergravity by compactification on a five-sphere. Its 
AdS§ solution is dual, in the planar limit and at strong 't Hooft coupling, to d = 4, M = 4 
super Yang- Mills (SYM) theory, which is a conformal quantum field theory Q . 

The study of fluctuations around other bulk configurations, which are interpreted as 
the duals of Renormalization Group (RG) flows of M = 4 SYM theory driven by relevant 
operators, is interesting, because it yields the correlation functions of the respective dual 
quantum field theories in a regime that is not accessible by ordinary perturbation theory. 
For example, the GPPZ flow which is a super symmetric mass deformation of Af = 4 
SYM theory, possesses many qualitative features of pure N = 1 SYM theory, in particular 
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quark confinement. The two-point functions in the GPPZ flow, which have been calculated 
using holography, exhibit a spectrum of particles, which must arise from as yet unknown 
non-perturbative effects in the field theory. These particles are interpreted as glueballs of 
the M = 1 SYM theory, which lies at the infra-red end of the RG flow ||]. Naturally, a 
calculation of the glueball scattering amplitudes would be very desirable. 

Until recently, the holographic calculation of n-point functions with n > 2 involving 
the operators that are dual to the active scalars 2 of holographic RG flow backgrounds 
was considered unfeasible, because their fluctuations couple to the fluctuations of the bulk 
metric even at the linearized level. In contrast, inert scalars do not couple to the bulk 
metric at the linearized level, and the calculation of the corresponding three-point functions 
is rather straightforward 0. In a recent paper ||], however, enormous progress was made 
exploiting a gauge invariant formalism, in which the true degrees of freedom of the bulk 
metric decouple from the active scalar fluctuations at the linearized level. Thus, this 
formalism simplifies the calculations of two-point functions for the active scalars and makes 
higher n-point functions accessible by the use of the Green's function method. This was 
demonstrated for the case of a single active scalar, which is dual to an operator O, by 
calculating the three-point function (OOO). An important result of ^ was the proof 
of Bose symmetry, which was not obvious from the formal expression resulting from the 
Green's function method. This proof involved two steps, which we shall summarize here. 
In the first step a field redefinition was used in order to eliminate those terms from the 
quadratic source of the field equations that contain two radial derivatives. Then, after 
using momentum conservation, some terms in the radial integral were integrated by parts, 
and the final, Bose symmetric result was obtained. Moreover, it was shown (for the GPPZ 
flow) that the boundary terms from the integration by parts cancel the contribution from 
the field redefinition. 

In this paper, we continue the programme of [|| considering a system of one active 
and one inert scalar coupled to bulk gravity in a generic holographic RG flow background 
and finding the expressions for all non-local three-point functions of the respective dual 
operators. Applied to the GPPZ flow, the results are then used to calculate the scattering 
amplitudes for the glueballs generated by these operators. Thus, for the first time, we are 
able to predict glueball scattering amplitudes for an M = 1 SYM theory using holography. 3 

In order to avoid being repetitive we chose not to draft this paper in a self-contained 
fashion, but to make essential use of the material presented in |8|. Hence, we urge the reader 
to consult that paper first in order to become familiar with the gauge invariant method 
and our notation. We also refer to j|] for more references to the relevant literature. 

An essential feature of pure Af = 1 SYM theory is missing in the GPPZ flow, namely the gaugino 
condensate. There exists a family of analytic solutions of bulk backgrounds, which include a non-zero gaug- 
ino condensate [M, but the linearized fluctuation equations around these backgrounds are not analytically 
solvable, so that one must resort to numerical methods to obtain, e.g., the mass spectrum of states H. 

2 In the common nomenclature, an active scalar is dual to the operator driving the RG flow and has a 
non-zero background value, whereas the other scalars are called inert. 

The qualitative behaviour of glueball scattering in confining gauge theories with string duals has been 
discussed before, see the recent talk M and references therein. 
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The rest of this paper is organized as follows. In Sec. || we review the equations of 
motion using the gauge invariant formalism developed in ||. Sec. [3|, which contains the 
main achievements of our work, is devoted to the study of three-point functions in a generic 
holographic RG flow background containing one active and one inert scalar field. Together 
with the bulk graviton, the three respective dual boundary operators give rise to a total 
of ten independent three-point functions, all of which will be calculated. We proceed as in 
Q performing first a field redefinition and then a suitable integration by parts in order to 
render the final results Bose symmetric. However, we shall not be concerned about whether 
the boundary terms resulting from the two steps cancel each other, as the net boundary 
terms would contribute only contact terms to the three-point functions. Our final expres- 
sions for the three-point functions have the form of generically divergent integrals. The 
divergences, which have the form of contact terms, can be easily understood and removed 
on a case-by-case basis by comparison with the results of holographic renormalization, 
and, thus, we shall not do it explicitly for the general setup. Such an attitude is justified, 
because contact terms have no effect on the scattering amplitudes. 

As an application of our results of Sec. |3] we shall analyze the GPPZ flow in Sees. [I| 
and [|. In Sec. ||, we start by repeating the analysis of the bulk-to-boundary propagators, 
which give rise to the two-point functions. Then, we shall present the expressions that 
encode the non-local three-point functions and scattering amplitudes. Moreover, we shall 
illustrate how the divergences of the three-point function integrals are understood and 
removed by using holographic renormalization. Finally, in Sec. |f| we discuss some physical 
interpretations of the calculated scattering amplitudes, in particular the possible glueball 
decay channels. Some useful relations for the GPPZ flow can be found in appendix |A], and 
numerical results for the decay amplitudes are listed in appendix [B[ 



Let us stress that we are interested only in the non-local three-point functions. A 
detailed analysis of the contact terms would be a worthwhile exercise, because they play a 
crucial role in the consistency of the subtraction procedure [|l(]] and contribute to certain 
sum rules [11|. Such an analysis would have to take into account not only the contact 



terms that result from holographic renormalization [12, 10, [L3|], but also the use of the 



gauge invariant fields, the field redefinition and the integration by parts, which occur 
in our calculation. In particular, analyzing the contributions from the use of the gauge 
invariant fields might be rather challenging, because one would need to translate the fields 
into the axial gauge, in which holographic renormalization has been carried out. Moreover, 
this translation must be done to third order in the fluctuations (to obtain all local terms 
in the three-point functions), but we have only derived the linear relations, which are 
sufficient for the analysis of the bulk equations of motion. Alternatively, one might try 
to perform holographic renormalization in a gauge invariant fashion. However, contact 
terms can be calculated only on a case-by-case basis, e.g., for the GPPZ flow. Thus, it 
would be impossible to include them into our general results of Sec. [|. Having in mind to 
calculate the scattering amplitudes for on-shell glueballs in the GPPZ flow, we feel justified 
in omitting them completely. 
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2. Field Equations 



In this section, we shall present the equations of motion to quadratic order in the fluctua- 
tions for a bulk system containing one active scalar, (f>, and one inert scalar, a, coupled to 
gravity. Our main reason for this restriction of the number of scalars is that, generically, it 
is impossible to diagonalize the effective mass terms of scalars of the same kind (active or 
inert) in holographic RG flow backgrounds. Thus, we ensure that the linear second order 
ODEs for the true degrees of freedom are not coupled, which makes the use of the Green's 
function method straightforward. Our results are easily generalized to the case of several 
inert scalars with diagonal effective mass terms. 

The equations of motion shall be presented in the gauge invariant approach of |8| . In 
this approach, the field fluctuations are combined into a set of gauge invariant variables, in 
terms of which the equations of motion are expressed. The remaining fluctuations describe 
gauge artifacts and are explicitly dropped. 



For completeness, we shall review in subsection 2.1 the general relations defining a 



holographic RG flow background and the definitions of the gauge invariant variables. Then, 



in subsections 2.2, 2.3 and 2.4, respectively, the equations of motion for the active scalar, 
the inert scalar and the traceless transversal parts of the bulk metric fluctuations are 
presented. 

2.1 Preliminaries 

A generic holographic RG flow background with one active scalar satisfies [14, |D|, [16| 

ds 2 = dr 2 + e 2A{r) i] ij dx i dx j , 
d r A(r) = -J^Wtf) , (2.1) 
d r 4 = W^4>) , a = 0, 

where the superpotential TV is a function of the active scalar and satisfies the following 
functional equation, which is imposed by the supersymmetry of the background, 

\ W l ~ d~L W2 = Vla=0 ' (2 ' 2) 
Here and henceforth, derivatives with respect to the fields are denoted by subscripts, as 
in Wfj, = dW/d(f). Notice that, in general, fl2.2j ) does not include the inert scalar, although 
in some cases {e.g., for the GPPZ flow) it might be possible that the full potential V is 
expressible in terms of a suitable superpotential containing also the inert scalar. 

The couplings between the inert scalar and the active scalar are strongly restricted by 
the fact that the relation V a = dV/ da = holds in the background due to the equation of 
motion for a. Taking a derivative of this relation with respect to r we find 

= d r V a = V^W* , (2.3) 

so that we find also V a $ = 0. Applying this argument recursively one arrives at the 
conclusion that 

= (2.4) 
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for any number of (^-derivatives. Notice that this argument does not hold in the more 
general case with several active scalars. In that case the right hand side of ( |2,3[ ) contains 
a sum over all active scalars, and the members of the sum need not be zero individually. 

We shall now review the definition of the gauge invariant variables. The bulk system 
is treated in the time slicing formalism, where the bulk metric is written in the form 

ds 2 = (n 2 + riin^dr 2 + 2n i drdx i + g ij dx i dx j . (2.5) 

The fluctuations of the various fields around the background (^T]) are introduced by 

4> = 4>(r) +ip , a = a, 

rii = Vi , n = 1 + v , (2.6) 

9ij = e 2A(r) (r]ij + hij) . 
Furthermore, the metric fluctuations /i*- are split as follows, 

h) = h TT ) + d% + V + + _La}fc , (2.7) 

where /i TT * denotes the traceless transversal part, and e l is a transversal vector {die 1 = 0). 

The inert scalar, u, is gauge invariant to lowest order. Therefore, we shall use the same 
symbol to denote the corresponding (all order) gauge invariant variable. The remaining 
first-order gauge invariant combinations of the fluctuations are 

a = ^ + w <t>W' (2 - 8) 

h = » +d \w)> (2 - 9) 

c = d^ + nJ^-^e 2A d r H , (2.10) 

cf = Uy - <? A d T i , (2.11) 

e} = h TT ) = Up{ . (2.12) 

Here and henceforth, II*- and IK? denote the transversal and the traceless transversal pro- 
jectors, respectively, 

& f)- 

n* = 5) - , (2.13) 

1 / \ 1 

jik _ ± /TT*fcrr .. i rrmfcl rrirrfc 



a)* = - (n*^, + qn* j - j—^)^ • (2.14) 

The gauge invariant approach is embodied in a simple recipe, according to which one 
expands the field equations to the desired order and replaces the fluctuations as follows, 

Cxi 

ip -» a , v -> b , v l -» cf + — , ^ e} . (2.15) 

The fluctuations e l , h and if represent gauge artifacts and are explicitly dropped, while 
the inert scalar a remains unchanged. 

Finally we mention that, as in Q, we raise and lower the indices of fluctuations and 
of partial derivatives using the flat metric. 
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2.2 The Equation for the Active Scalar 

The equations of motion for the active scalar can be taken over from Sec. 4 of [|| with only 
minor modifications, which come in the form of additional terms that are quadratic in the 
inert scalar. Terms with only one a cannot occur because of (|2.4|), For completeness, we 
repeat the full expressions here. First, the scalar equation remains 

2d 



W3 r + e~ 2A n — Vfhd) I a — W^~' ZA c - WM - 2VJ> = J a 



2A, 



where the source J a is now given by 

Ja = \y<t»t><t>a 2 + V^b 2 + 2V H ab - W^b + (d r a){d r b) + ^W^d^ 



+ e~ 2A 
- W s I < l 



■2hDa - (0'h)(0,u) + <■<),.„ r 2 ( ,t + — ) <),() r o 



□ 



dd 3 

+ eUd i cPa-W lt> d l d J -W 6 1 



, \ v 2 

"3 \ ' ' " <P q "J 2 ^" Ja ' 

Second, the normal component of Einstein's equation is 

-AWe~ 2A c + AW^a - AV^a - 8Vb = J c , 
where the source J c is easily generalized from (4.6) of ||, 
j c = AVb 2 + 8V±db + 2V H a 2 - 2{d r a) 2 + (e~ 2A c) 2 + 2<r 2A {d i a){d i a) 



+ AW 6 e 



-2 A 



□ 



did + 2We l i d r ej - AWe~ 2A e l A I d 3 + 



dh 

□ 



- jidre^idrej) + e~ 2A [ did? + ^-c ) d r e 



did 3 



□ 



— e 



-4 A 



-2 A 



8% 



-{dd'Wdj) + -(did^djd 1 ) + 2{did 3 ) D 



-c + 



1 

□ 



1 

□ 



ejde? + ^e£)(^) - ^ejX^e}) 



+ 2V aa a z - 2(d r a)(d r a) + 2 e - 2/1 (dV)(d^) . 
Third, the mixed components of Einstein's equation yield 



1 



e~ 2A adi - 2Wdib - 2W (t> d i a = J { , 



and the source </,■ now is 



Ji = -Wdib 2 + 2{d l a){d r a) + e- 2A (n^c)(^6) + ^{djb)(d r el 



2A d j c Uc j 
□ 1 



~ le- 2A (d je l)(&d k - d k d>) - \e~ 2A dp4 - ^ e - 2A (d,b)(d 3 d l + d t d 3 ) 

+ 2{dio){d r o) . 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 
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As in ||, we first define a = (W/W^a and solve (|2.18|) and (|2.20|) for b, c and dj, with 
the results 



1 a ■ 

b = %a-— , 

W 2 2W □ 

□di = -2e 2A Upj , 
W 2 



(2.22) 



e" 2A c 



Then, we substitute ( |2,22| ) into ( |2,16| ) , which yields the following second order ODE for a, 



D + e~ 2A n) a = J- a 



We have abbreviated 



( W l 



d 



and the source term in ( |2.23| ) is given by 



a ~w/ a 4 Jc 2 



W d-l 



W 2 



-w 



Or 



(2.23) 



(2.24) 



8 r + 2 \W 66 - 



d 



-w 



□ 



(2.25) 



W d-l 

In the source (2.25), the first order terms of the solutions ( |2.22| ) should be substituted. 
2.3 The Equation for the Inert Scalar 

The equation for the inert scalar is quite easily obtained from ( [2.16 ) and ( |2,17 ) by dropping 
the terms with and and generalizing the derivatives of V. One obtains 

2d 







Wd r + e- 2A D-V arr )a = J a 



d-l 



(2.26) 



where the source J„ is 



Ja = -^VaaaCT 2 + Vaa<f><?a + 2V aa cb + (d r Cx)(d r b) 



+ e" 



-2A 



- '2 b D n - (0'b)(0,a)+ r!),.a + 2 ( <? + — ] i),O r a 



(2.27) 



+ d r d l + d r — &<r-2 



'^4w(d i + ^)n i a + c',i) l i)>n 



□ / " d-l 
2.4 The Equation for the Graviton 

The tangential components of Einstein's equation provide the equations of motion for the 
traceless transversal modes e*-. Following the notation of ||, we first write these components 
in terms of time-slice hyper surface quantities, which gives 



- d r {nK)) + n k V k {nlC)) + {nK))(n1C k k + d r lnn- n k d k Inn) 

+ nV'djn + {nlCi)V j n k - {nK, k )V k rj + n 2 R) + 2n 2 S) = , (2.28) 
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where 



S)=g lk {d k ct>d^ + d k ad j a) + 



-5)V(cP,a) 



(2.29) 



d-1 r 

In order to obtain an equation containing only the gauge invariants e*-, a and <7, we should 
expand ( 2.28| ) to second order and substitute the solutions ( 2.22j ) for b, c and d\ However, 
it is clear that the trace and the divergence of the resulting equation would both vanish 
by virtue of the equation for a, ( [2.23 ), because the scalar equation is implied by Einstein's 
equation via the Bianchi identity. Thus, we can project onto the traceless transversal 
components using the projector ( 2. 14] ) in order to find the only independent equation that 
is still missing. The result is 







2d 
d-1 



Wd r + e' 2A a 



J 



j > 



where the quadratic source terms are given by 



(2.30) 



jj = Uf l {(d r e l m )(d r ef)-e- 2A 



e™{2d l d m e n k - d m d n e l k ) 



+ ^d l e™)(d k e n m ) + (d m e' n )(d n ef) - (d m e n k )(d m e l n ) 



I \ I an „m N i 



w 

+ 2|'3t 

1 w 



+ 



Or 



w 



□ 



{d m d r e l k ) ( — d r d ) + e- 2A da e { 



□ 



■d r d 



d r a ) + e~ 2A d l dd k d 



□ 



d l d k 

i 

□ 



d r d — 4e' 



-2 A 



w 



w 



ad r e{ (2.31) 



d l dd k a + d l ad k a 



3. Three-Point Functions 



3.1 General Form 

The dynamics of the bulk fields, which is governed by the equations presented in the 
previous section, encodes the two- and three-point functions of the dual operators. For 
completeness and in order to outline our conventions for the presentation of the final 
results, we shall in the following review how the correlations functions are obtained from 
the (sub-leading) asymptotic behaviour of the bulk fields. 

The gauge invariant fields a, a and e*-, which are the true degrees of freedom of our 
bulk system, are dual to the operators O, S and Tj, respectively. Here and henceforth, we 
denote by Tj the traceless transversal part of the boundary energy-momentum tensor, 



(3.1) 



The duality between the bulk fields and the boundary operators is made explicit through 
the couplings 



d d x 



a(x)0(x) + a(x)E(x) + -e){x)Ti(x) 



(3.2) 
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where a, a and e* are the prescribed asymptotic boundary values of the bulk fields a, a 
and e % p respectively. 

The other components of the energy momentum tensor are not independent, because 
the (anomalous) Ward identities Jl^, [l3| imply the operator identities 

di@) = , Q = (30 , (3.3) 

which are valid in all correlation functions with distinct insertion points. For vev flows we 
have (3 = 0, whereas in operator flows (3 = —{d — A)</>, where A is the conformal dimension 
of O at the ultraviolet conformal fixed point of the RG flow. 

In order to start, let us denote by tp a generic bulk field. Its behaviour in the asymp- 
totically AdS region of the bulk space-time is described by the generic expansion 4 

iP(p, x) = A )/ 2 $( x ) + • - - ] + p A/2 $(z) + •••]. (3.4) 

Throughout this section, we shall use the variable p = e~ 2r , so that the asymptotic region 
is given by p — > 0. Here, A denotes the conformal dimension of the operator dual to 
tp, which we shall call Vl/. 5 The ellipses in ( |3.4D stand for the sub-leading terms in the 
two series, which are power series in p, whose coefficients depend locally on tp and tp, 
respectively. If A — d/2 is an integer, the leading series contains also logarithms. The 
function tp is called the response function, and its non-trivial dependence on the source tp 
stems from the condition of regularity of tp in the bulk interior. Moreover, the response 
function tp determines the exact one-point function of the operator ^. More precisely, we 
have 0, 0, 0] 

(*) = (2A - d)tp + contact terms , (3.5) 

where the contact terms are finite, but in principle scheme dependent. We shall not be 
concerned with the contact terms in this paper. A slightly different version of (^^) is 
needed for the exact one-point function which is 

<^> = ~(2A-d)4 = ^. (3.6) 

Here, contact terms have been omitted. The factor 1/4 stems from the 1/4 in front of the 
Einstein-Hilbert term in the bulk action in our conventions § . 

Let us focus our attention now on the response function, tp. In order to determine it 
we first observe that tp satisfies a second order ODE of the form 

(V 2 - M 2 ) tp = , (3.7) 



4 In the general analysis, we shall set the length scale of the asymptotic AdS region to one for simplicity, 
i.e., L — 1. 

5 We restrict A by d/2 < A < d. The upper bound means that we consider relevant or marginal 
deformations, whereas the lower bound stems from the standard AdS/CFT prescription. An extension 
down to the Breitenlohner-Freedman bound p^ , i.e., for d/2 — 1 < A < d/2, has been developed in 
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where M 2 is an effective mass term, V denotes the background covariant derivative, and 
J^p is a higher order source term. Thus, after defining a covariant Green's function by 



V 2 -M 2 )G(z,z') = ^=^, (3.« 

V9(z) 



the general solution of Q3.71) is given by 



${z) = / d d yK(z,y)iP(y) + / d d+1 z' y/§^)G(z, z')J^(z') . (3.9) 



Here, z is a short notation for the variables (p,x), and x and y are boundary coordinates. 
Notice that the bulk integral is cut off at p' = e, and also that p > e, because of the 
regularization procedure. Moreover, K(z, y) denotes the bulk-to-boundary propagator of 
the field ip. 

We are interested in the near-boundary behaviour of ip, and it is very helpful that in 
asymptotically AdS spaces the Green's function asymptotically behaves as 

G(z,z')tt-J^ 1 K(x,z') + --- . (3.10) 

Hence, let us consider the field at the cut-off boundary by setting p = e, and let us 
switch to momentum space, where we can use momentum conservation of the propagators, 
K(p,p; q) = K p (p)5(p + q). Then, (gj) leads to 

iP(e,p) « K p (e)#(p) - J ^ dA ^K p {p)J^p,p) • (3.11) 

From this general expression one must extract the response function ip, which represents, 
by virtue of ( |3.5[) , the exact one-point function of the dual operator The two-point 
function can be read off easily from the asymptotic behaviour of the bulk-to-boundary 
propagator K p . 

The various three-point functions that contain ^ are obtained from the interaction 
integral in ( 3.1 1| ) . In order to be more explicit, let us consider a generic three-point func- 



tion (\E f i^ r 2 1 3 / 3) 5 where the n = 1,2,3, can be identical or different. Their respective 
dependence on the momenta p n is implied. This three-point function can be obtained from 
the integral in ( 3.1 1| ) setting, e.g., ip = ifii. The term in that is responsible for the 



three-point functions (^ / i^ , 2^ / 3) has the generic form 

Jil>l(P,Pl) = J dP2dp35(p 1 +p2+P3)X(pi,-p2,-P3)K 2 (p)K 3 (p)ip 2 (-p2)'4>3(-P3) , (3.12) 

where we have substituted the linear solutions for the bulk fields ip2 and ip^, and X is an 
operator containing also derivatives with respect to p (or, equivalently, r) acting on the 
bulk-to-boundary propagators K 2 and K%. It is also important to notice the minus signs 
in front of P2 and ^3. Differentiating ( 3.12| ) with respect to the sources ip2 and 1P3 and 



substituting the result into ( p. 11 ) yields the three-point function of the form 

(*!* 2 * 3 ) = -5( Pl +p 2 +P3) j dre dA X^K X K 2 K Z , (3.13) 
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where X123 is again an operator acting on the bulk-to-boundary propagators, which can be 
obtained easily from the source . 

Eqn. ( |3.13D is the general expression for the three-point functions, and our main work 
will consist in calculating the operators ^123- To explain this formula more precisely several 
clarifications are in order. First, the factor (W^/W) 2 has been inserted for continuity with 
our previous calculations in [||], where it automatically appeared when using the field a 
and its corresponding bulk-to-boundary propagator, K. In fact, in the case = O, the 
operator A123 can be read off directly from the source J a . In the other cases, ^1 = X and 
^1 = Tj, we must multiply a factor iW^/W) -2 to the operators read off from J a and Jj, 
respectively, in order to compensate for the factor in ( 3.13| ), Second, we must correctly 



carry out the functional derivatives with respect to the sources, if \E , 2 and ^3 are identical, 
which symmetrizes the operator of interest with respect to the indices 2 and 3. Third, 
because of ( p\^ ) and (|3.6|), we must include a factor 1/4, if ^1 = T' 1 -, and a factor 2 for 
each at the positions 2 and 3. Forth, the reader is reminded of the minus signs in front 
of P2 and P3 in ( 3,12j ). Last, the integral in ( 3.13| ) is divergent in most cases. This is not 



much cause of concern, because the divergences can be predicted and subtracted using the 
counterterms of holographic renormalization. This is done best on a case-by-case basis, 



and we shall illustrate it for one of the three-point functions in the GPPZ flow in Sec. 4.3 



3.2 The Operators ^123 

We shall, in this subsection, calculate the operators X123 for all ten independent three-point 
functions of O, X and T % -. For this purpose, we take the equation of motion for the dual 
field of one of the operators involved and consider the source term containing the other 
two fields. If the source contains terms with radial derivatives on both fields, we use a 
field redefinition to eliminate those terms. Then, as explained in the previous subsection, 
the operator X123 can be read off easily from the new source. It will contain terms with 
radial derivatives (with respect to r) acting on one of the bulk-to-boundary propagators 



in the integral ( 3.13 ). These derivatives are denoted by d\, 82 and 83, which act on Ki, 
K2 and K3, respectively. In most cases we need to integrate by parts some of the terms 
in order to obtain a Bose symmetric expression. The resulting boundary terms (as well as 
those from the field redefinition) can be dropped, because we are not concerned about the 
scheme dependent local terms of the three-point functions. 

The three-point functions involving operators of different kinds can obviously be cal- 
culated in more than one way. We have actually performed all 18 possible calculations and 
cross-checked the results for the mixed three-point functions. In the following, we shall 
present the ten independent expressions starting with Xooo, which we repeat from for 
completeness. 

The Correlator (OOO) 

This three-point function can be calculated only by considering those terms of the source 
J a of ( 2.23] ) that are quadratic in a. These are obtained by inserting the linear solutions 
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for b and c into the general expression ( 2.25Q dropping d{, e*- and a. Thus, we find 



Ja 



1 fw d 



4 V W 



2 S (afnia!)-^ 919 ^' 



□ 



□ 



-a 



+ 2 



W^\ 2 &a! 



W / □ 



did,' 



1 

+ 2 



<9r 



re 



□ 



are a' + 



9 



9* a' d l 

did, (a! dia} — ac 



1 fw d 



2 v re 



re,, re ( 



□ □ 
2" 



re re 2 



a 



2an a + ^ [(5f5)(na)] - ^(^S)(^a) 



(3.14) 



where D 2 is the second order differential operator defined in ( p. 24 ) , and we have abbreviated 



a' = d r a. 



As in [Q, in order to facilitate the integrations by parts, we remove the terms of J a 
with r-derivatives on both fields by the field redefinition 



a — ► a + 



1 /re. 



8 V re 



2^- [allla 



8 i 8- 



+ 



re 



d l a 
□ 



(3.15) 



In terms of this new field, the source becomes 



d i d j ( ,didi 

a a 

a \ □ 



□ 



□ 



Or 



re, 



re 



d l a! d l a . <9* , . N , 

c^a + 2-^-aja + — (a c^a] + a c 



□ 



□ 



1 /re, 



2 v re 



-2A 



ana _ ^1 [adiCPci) + ^{dfa){d l ~a) 



+ ( ^)V 2A |2aDa + I ma)(na)] - -(^(^a) - 2 f^fi] (/W) 



1 

+ 2 



£ 2 



re re 2 



a 2 



D 1 



W± 

re 



re^ 
re 



a 2 - 2- 



• 3 1 a -a I + 



□ 



□ 



did 

□ 



l d,; _\ ( d^d., 



□ 



□ 



(3.16) 



We can read off X123 from ( 3.16| ), where we have to take care to symmetrize the 
indices 2 and 3 because of the functional derivation with respect to the sources. Using also 
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momentum conservation, pi+p 2 +P3 = 0, we find the operator 

4 



Xooo 



1 



Or 



w 

(P2 • P3? 



03 H — 9 — C2 H — — °1 



PlP 2 



P2P3 



+ 



r 



w 



1 



1 



2 2 z 1 2 2 

P2P3 



1 



P2 n 



P2 ■ P3 -o + -9 #L + Pi • P3 -o + -9 # 



1 



pi p§ 



"Pi • P2 -5 + -2 ) #3 - P2 • P3 ( -9 + "a (#1 + + 93) 

\Pi pi J \p 2 pi 

4 



+ 



+ 



1 /w* 
w 



-2 A 



\{pI+pI+pI)+ . m 



(Pi -P2) , (Pi -P3) , (P2 -P3) 



+ 



+ (Pl -P2) 2 



-2A 



1 



(P2 • Pzf 



1 



1 



o2 ^2 



+ (Pi -P3) 2 



1 

— 2 2 
Pi P2 



1 

+ 4 



+ <^ 2 



W± 
W 



w 



vP2 P3 

2(p?+pl+p§) 

(P1-P2) 2 , (P1-P3) 2 (P2-P3) 2 



2 2 
P1P3 

1 t 1 

— 2 2 
Pi P3 



+ 



P2P3 



p\p\ 



+ 



P\P\ 



p\pI 



w^V 1 fw^ 1 



w 



w 







]}- 









P2 -P3 



1 1 

— 2 2 
P2 P3 



(3.17) 



Because of the identity 



(di + d 2 + 83)^X2X3) = d r {K x K 2 K3 



(3.18) 



we can integrate these two terms of ( |3.17 ) by parts in the integral in ( 3,13j ) taking note 
also of the identity 



w 1 " 



w 



D 2 . 



(3.19) 



Hence, the last term in ( 3.17| ) is cancelled, and the minus sign of the last term on the 
penultimate line is reversed, rendering the final result totally symmetric in the indices 1, 2 
and 3, as it should be. Thus, the final result is 



Xooo = 2 



Or 



w 



(Pi -P2 f 
P\P\ 



03 + 



w 



Pi • P2 ( ~2 + \ ) 9 3 

.Pi Pi, 



1 (W< 



+ 



w 



-2 A 



D 1 



w 



1 2 , (Pi -P2) 

0P1 + Tl— 

2 p{p l 2 

(Pi -P2) 2 1 ' 

p\p\ 3 



+ 



D 



w 
w 



-2 A 



(pi • P2) — + — 
Pi p 2 



1 



2p? 



w*\ 2 1 rw, x ! 



w 



w 



+ cyclic . 
(3.20) 
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The Correlator (OOT) 

In order to calculate the correlator (OOT), we consider the terms of J a that are of the 



form ae l y These are 



Ja 



1 fw^ 
w 



did 3 f~, v 
□ V 3 



did 3 _, 



□ 



a' e • + 



□ 



□ / 3 



1 

+ 4 



'■' W 

2 



□ 



□ I 3 



1 (W± 
W 



e -2Adidi ,~ d js + e -2AAg. dj ~ 



(3.21) 



The terms with two r-derivatives are removed by the field redefinition 



a 



a+ 8[w 



did 3 
□ 



did 3 _\ i,di 

a e'- H 

a 3 □ 



—a ) De!- 



(3.22) 



which leads to the new source 

2 



Ja 



1 



d 



did 3 



did 3 / ,/ . 



□ 

1 (W± 

w 



-2A 



— [d k ad ej 



dd 3 



□ 



1 



D 1 



W0 

w 



did 3 
□ 



(5e$) 



*/ 

a e,- — 3 — 

a 3 □ 



—a De 
□ 



□ 



d, 



d k e) + ^ 



a et- H 

□ / 3 □ 



<9 J 

—a I Dell- 
□ 



+ e" :L4 e!-d i <9 J a 



(3.23) 



From ( 3.23| ) it is straightforward to read off the operator Xoot, where we must remember 
to multiply by a factor of 2 for the T insertion. Using momentum conservation, it can be 
written in the form 



-2A 2 2 
PlP2 



+ 



dr 



W 



[(pl-pl+pl)(di + d 2 + d 3 ) 



+(£ -Pi- Pz)di + (p| -pj- Pl)d 2 + 2(pl -v\- p 2 2 )d 3 ] 



+ 



W 



(pI-Pi+p!) 



(3.24) 
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Integrating the term containing (81+82+83) by parts we cancel the term on the penultimate 
line and end up with 



lyr ikPlkP 2 L-aA. 



V l TT IKi L^r< 



. 2 2 
P1P2 



Wa 



M -Pl- Pl)8x + (vl -Pi- Pl)8 2 + 2(pl -pj- p 2 2 )8 3 



! ( ] e~' A (pt +p$+pt- 2pfpi - 2pfpi - 2p z 2 pi) 



Obviously, the operator is symmetric in the indices 1 and 2, as it should be. 



(3.25) 



The Correlator (OTT) 

For this correlator we need to keep the terms in J a that are quadratic in e*-. It is straight- 
forward to find 



w 



+ 4 e " 



8 l 8 



2A I u U J Jr-n^k 



□ JP<) + ~^<W^) - -(d ie {)(8«e}) + -- l(8 t e{)ne« 



The first term is removed by the field redefinition 



1 



a -> a + ~n}(elef) , 



(3.26) 



(3.27) 



which leads to the new source 

1 



Ja 



1 



-2 A 



_> _ (de J e k 



1 



U8 l e{We k j ) + "- (d ie ? k )ne 



8'' 



(3.28) 



From ( 3.28 ) we easily read off the operator (remember the factor 4 for the two T insertions) 

W l\ 

- -r^- I (pi • P2&3 + Pl • P38 2 ) 

(3.29) 



v ik TT iniTT kn 

^OTTji — ~2 LL 2 jn LL 3lm 
Pl 



1 



+ ^~ 2A (P\ +Pt +Pl ~ 2p\p\ ~ 2p\p\ - 2p\p\) 



This operator is symmetric in the indices 2 and 3 because of the functional derivation with 
respect to the sources S. 
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The Correlator (COS) 



Due to (2.4), there is no term of the form aa in the source J a- Hence, this correlation 
function vanishes, 

(COS) = . (3.30) 

The Correlator (OSS) 

The terms of Ja that are quadratic in a are 

/ Wl \ d i 1 / W \ 

J a = -2 \W H -^J- {v'dicx) + - [ W V<t>«« ~ v 2 

+ e- 2A l [ ^[(d l a)(na)]-^a)(d l a) 
As terms with two r-derivatives are absent, we can directly read off the operator Aoee, 



W 2 

Wtj, v p{ 

1 



W l\ 

- — Upi ■ p 2 d 3 + pi ■ p 3 d 2 ) 

' ' (3.32) 



+^~ 2A (pi +pl+pl- 2p\p\ - 2p\p\ - 2p\p\) 
The Correlator {OUT) 

As there is no term of the form e*cr in the source J a , this correlation function vanishes, 

(osr;) = o . (3.33) 

The Correlator (EES) 

This three-point function can only be obtained by considering those terms of the source 
J rj of ( 2.26j ) that are quadratic in a. There is only one such term, which is 



Ja = ^V aaa a 2 , (3.34) 



from which we directly read off 



= (jf) 2 V m . (3.35) 



Remember that the factor (W^/W) is needed for the general convention ( 3.13 ). 
The Correlator (SET) 

For this three-point function we consider the term in J a that is bi-linear in e l - and a, which 
is 

J rj = e- 2A e)did j a , (3.36) 
from which we can easily read off the operator (remember the factor 2 for the T insertion) 

X^t) = 2 (^f) 2 e- 2A n 3 ;.W 2 • (3-37) 
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The Correlator (ETT) 

As there is no (e* ) 2 term in J a , this correlator vanishes, 



. 



(3.38) 



The Correlator (TTT) 

This three-point function is obtained from those terms of the source Jj of ( 2.3C| ) that are 
quadratic in rf. These are 



J) = U$ (eL)'(en'-e 



-2 A 



e™(2d l d m e n k -d m d n e 



n„l 



+ 2^0(^0 + (d m e l n )(d n eT) - (d m e l n )(d m e n k ) 



After the field redefinition 



we obtain the new source 



Jj = -n j7 e 



e j + 2 i l m k 



e™{2d l d m e n k - d m d n e[) + -(d l e™)(d k e n m ) + (d m e l n )(d n e^) 



(3.39) 



(3.40) 



(3.41) 



Remembering to include the factors (W^/W)^ 2 for our convention, 1/4 for the first T and 
4 = 2-2 for the other two Ts, the operator Xttt is easily found from ( |3.41| ) , 

-2 



V ikm 
<*TTTjl n 



-2 A 



w 



nil' tt kk 1 tt mm 1 
ljj' LL 2W ii 3nn' 



X 2 (pi k >P2 5$' 6 l m , + Pim'PzSgSv + P2 m iP3${< S ii') 

Obviously, this operator is Bose symmetric. 



(3.42) 



4. GPPZ Flow Correlation Functions 

In this section, we shall apply our general results to the GPPZ flow with the aim of calcu- 



lating the glueball scattering amplitudes. To begin, we shall repeat in subsection |4.1| the 
calculation of the bulk-to-boundary propagators, which also encode the two-point func- 
tions, as they will be needed for the explicit calculation of scattering amplitudes later. The 
mass spectra and decay constants of the associated glueball states are read off in a stan- 
dard fashion. Then, by isolating the on-shell poles of the bulk-to-boundary propagators 



the external legs of the Feynman diagrams are explicitly amputated. In subsection 4.2, we 
provide the expressions for the non-zero three-point functions in the GPPZ flow. These 
automatically encode the three-particle scattering amplitudes by amputating the external 
legs. It is a fortunate fact that the on-shell (amputated) bulk-to-boundary propagators are 
polynomials. Thus, the radial integral in the scattering amplitudes is elementary. Through- 
out this section, we use the radial variable u defined by u = 1 — p = 1 — e~ 2r , as is customary 
for the GPPZ flow. The asymptotic AdS region of the bulk space-time is given by u — > 1. 
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4.1 Bulk-to-Boundary Propagators and Two- Point Functions 
Active Scalar 

Let us begin with the bulk-to-boundary propagator for the active scalar. The linear equa- 
tion of motion ( 2.23 ) for a becomes (in momentum space) 

„2" 



u(l - u)d z u + (2 - 2u)d u 



P 



a = 



(4.1) 



This is a hypergeometric equation, whose solution, which is regular for u = 0, is readily 
found. We have to be somewhat careful with the normalization, because the bulk-to- 
boundary propagator for a, K = (W^/W)!^ , should satisfy the generic asymptotic ex- 
pansion with a unit coefficient of the leading term. Hence, we find the bulk-to-boundary 
propagator K 



KM - f r (1±£) r (^ I f 



1 + a 1 — a 



;2;u 



(4.2) 



with a = \J\ — p 2 . Its asymptotic behaviour is |2C] 



k p (u) 



V3 



1 + ^(1- u) ln(l -u) + -H°(p)(l -u) + 



The function H°(p), which is related to the two-point function by 

(0(p)0(q))=6(p + q)H°(p) , 

is given by 



H°(p) 



P 



3 + a 



3 — a 



^(2)-V(l) 



(4.3) 



(4.4) 



(4.5) 



The spectrum of poles becomes clear after rewriting H in a series representation using 

1 



the formula [21] 



f^ \y + k x + kj 



This yields 



y + k x + k / 
2k + 1 



2 ^kik + l^Akik + ^+p 2 } ' 
Thus, we find particles with the masses 

m 2 = 4k(k + l) , k = 1,2,3,... . 

The residues at the poles, which represent the decay constants [ p2|| , are 

\f°\ 2 = 8k(k + l)(2k + l) . 



(4.6) 



(4.7) 



(4- 



(4.9) 
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For an on-shell momentum, p 2 = —Ak(k + 1), we have a = 2k + 1, and the hypergeo- 
metric function in ( f4,2|) truncates to a polynomial, whereas one of the T-functions in the 
normalization factor has a pole. As in 1| we identically re-write ( |4,2| ) as 



a 



k p (u) = ^-r 



and amputate the external leg as 



2 j v ; du V 2 ' 2 ' ' 



(4.10) 



I f°\ 

K (v) — 

pW ~ p 2 + 4k(k + 1) 



3(2fc + 1)(fc+ ^(l-n)P^(2n-l) 



2A; 



+ regular , (4.11) 



where P^*'^ denotes a Jacobi polynomial of degree k [21], 



ra=0 V ' V ' 



(4.12) 



Inert Scalar 

Next, let us consider the inert scalar. The linear equation of motion ( 2.26j ) for a becomes 



u(l - u)d 2 u + (2 - 2u)d u - V — + 2 



(7 



o . 



(4.13) 



Again, we have a hypergeometric equation, whose solution, regular for u = and normal- 
ized properly, is 



k;{u) = r 



3 + /3\ Z3-/3 



Vl -mF 



1 + /3 1-/? 
2 ' 2 



;2;u 



with (3 = ^9 — p 2 . Its asymptotic behaviour is very similar to §L 



KZ{u) « 



1 + ^^(1 - u) ln(l - u) + ^ S (p)(l - u) + ■ 



where H^(p), which embodies the two-point function (££), is given by 



Using (|4.6|) , we re- write (|4.16| ) as 

(p 2 -8) 2 * 



H^(p) 



E 



2k + 1 



2 ^ *;(£; + l)[4(A;-i)(& + 2)+ p 2 ] 
Thus, we obtain the mass spectrum 

m | = 4(jfe - i)(jfe + 2) , fe = 1,2,3, .. . , 
and the respective residues are 

\ff\ 2 = 8k(k + l)(2k + l) . 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 
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Notice that the mass spectrum ( [4.18 ) includes a massless particle. 

On-shell, we have (5 = 2/c + l, and, as for the active scalar, the hyper geometric function 
in ( [4.14] ) truncates to a polynomial, whereas one of the T-functions in the normalization 
factor has a pole. Hence, after first rewriting ( |4.14 ) as 



kz{u) = r 



(4.20) 



we amputate the external leg as 



p 2 + 4(k- l)(k + 2) 



1) 



+ regular . 
(4.21) 



Graviton 



Last, we consider the traceless transversal part of the graviton, e*-. Its linear equation of 
motion ( p. 30 ) becomes (in momentum space) 



u(l - u)d 2 u + (2 - u)d u 



4 = 0. 



(4.22) 



Also this is a hypergeometric equation, and its regular and properly normalized solution is 



k;(u) = r 



4 + 7 



4-7 



with 7 = \J —p 2 . Its asymptotic behaviour is |2C] 
~ 2 pV + 4) 



iff(u)«l-^(l-u) 



1 



(4.23) 



1 - u) 2 ln(l - u) + -# T (;p)(l - u) 2 + 



4 v ' 32 
The function H T (p), which embodies the two-point function by 



T^p)T l k ( q ))=Uf l 5(p + q)H T (p) , 



is given by 



H T (p) 



p 2 (p 2 + 4) 
16 



4 + 7 



+ V 



4-7 



^(3)-V(l) 



(4.24) 



(4.25) 



(4.26) 



After re- writing ( [4.26 ) as a series using (|4.6|), we obtain 



H T (p) 



p 2 (p 2 + 



k + 1 



^k(k + 2)[4(k + l) 2 +p 2 } ' 
Thus, we find particles with the masses 

ml = 4(Jfe + l) 2 , ft = 1,2, 3, 



/, . . . , 



(4.27) 



(4.28) 
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and the respective decay rates are 



\f i 



T|2 



8k(k + 2)(k + 1) 



(4.29) 



These states describe spin-two glueballs, but we postpone a discussion of the spin structure 
to Sec. ||, where it is needed for the calculation of the scattering amplitudes. 

On shell, we have 7 = 2{k + 1), and the hyper geometric function in ( 4.23| ) truncates to 
a polynomial, whereas one of the T-functions in the normalization factor has again a pole. 
As before, we first re- write ( |4.23 ) as 



K(u) 



2 + 7 



7 



2 + 7 2-7 



(4.30) 



and then amputate the external leg as 
\f k \ 



k;{u) 



p 2 + 4(fc + 1) 



2(fc + l)(fc + 2) 
k 



+ regular . (4.31) 



4.2 Three-point Functions and Scattering Amplitudes 



In this subsection, we shall re-write our general results for the three-point functions for the 
case of the GPPZ flow. In order to fix convention, we shall write the three-point functions 
in the form 

(*!* 2 * 3 ) = -5(pi +p 2 +p 3 ) J du yusK^Ks , (4.32) 

where the operators ^123 are related to the operators X123, which were presented in sec- 
tion EO, by 



3^123 



3 4A(u) 



2(1 



w 



X 



123 



1 



u 



X 



123 



(4.33) 



The symbols d n 



n 



1,2,3, that appear in the expressions for 3^123 are now derivatives 



with respect to u acting on the respective K n . The integral in ( [4.32 ) is divergent in all cases 
but one. The degree of divergence will be indicated, but the discussion of the divergences 
and how they are removed is postponed to subsection |4.3| . 

It follows from (4.32) that the irreducible scattering amplitudes take the form 



M 



123 



/ 



(4.34) 



where K n , n = 1,2,3, denote the amputated bulk-to-boundary propagators, which are 
given by the expressions in the brackets in (PH ), ( fL2lD and (|4~3l|) . Since the amputated 
bulk-to-boundary propagators are polynomials in u, the integral in ( |4.34j ) is in general 
elementary. 

For the GPPZ flow, there are only six independent non-zero three-point functions 
involving 0, S and Tj. In addition to (COS), (OTT,) and (TTS), which vanish in general, 
also (EES) vanishes in the GPPZ background, because of V aaa = 0. The operators ^123 
for the six non-trivial three-point functions are listed below. 
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The operator yooo was calculated already in ||, but we shall repeat it here for 
completeness. From ( |3.20| ) and ( [4.33| ) we obtain 



yooo 



32 
27 



_4n 2 fl - u)^K^d 3 - 3u 2 Pl ■ p 2 ( * + * ] 0, 



+ -u(l-u) 



+2u(3u - 2 



p\p\ 

l ( 2 , 2n , (Pi -P2) 3 

4 



n P2 



3 

+ 4 n 



(Pi • P2Y ( -j + S 

pi pi 



(pi+pi) 



(gl -P2) ; 



+ 1t 



2(1 - u) + ~ 



+ cyclic . 



(4.35) 



Since K ~ \/3/2 to leading order, the integral ( [4.32 ) for (OOO) is finite. 
For the correlator {(DOT) we find 

2 



y OOT } = U ^lP^P l 2 9p 2 p 2 



6- 



1 



2 2 
-Pl?>2 



- 4u 2 [(p 2 - p 2 - p 2 )^ + (p 2 - p 2 - p|)a 2 + 2(p 2 - pi - p 2 2 )d 3 ] 
+ <p\ +p 4 2 +pI- 2p\p\ - 1p\p\ - 2p\p\) 



(4.36) 



The first term in the braces leads to a logarithmic divergence in the integral in (|4.32| ) 
For the three-point function (OTT) we obtain 

?il - JagDsft^L^JL- {-Au[(pl -p\ -p 2 )d 3 + (p 2 2 -p\ - p 2 3 )d 2 ] 
+(Pl + Pi + P3 " 2p?P2 - ipXpI - 2p|pi)} • 



(4.37) 



Again, the first term in the bracket leads to a logarithmic divergence in (4.32). 
Next, for (OSS) we find 



M2-3^ 1 



(1 — u) 2 3p 2 1 — u 
+(Pi + Pa +P| - 2p?P2 " 2p?P3 - 2^!)} ■ 



{-4u[(p 2 - p 2 - p 2 )d 3 + (pi -pX~ Pl)d2 



(4.38) 



The first term leads to a logarithmic divergence in the integral in ( 4.32 ), because of the 
leading behaviour K a ~ y/\ — u. 



The expression for (EST) is rather simple, 

y^ T ) = u 3 f lPlkP l 2 - 



j ^^-^ (1 _ u )2 • 

Again, the integral in ( |4.32 ) is logarithmically divergent. 

Finally, the operator for the correlator (TTT) takes the form 



(4.39) 



-\ 1 ikm 
yTTTjln 



kk 1 t-i mm 



, 1 



U 



l jf ii2ll> ^3nn' 2 ^ _ s 2 

2 (pwPzStf&m- +Plm>p(Sk>4 +P2m>P3%$) 
-I' , _ Ji'sn'd' , - ._/«»' rf' " 



(4.40) 



+Plm'i>2 b{A>+Plk>PW> & m ' +PlM 6 k' 5 m' 

This operator leads to a linear divergence in the integral in ( |4.32| ) . 
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4.3 Finding Divergences from Holographic Renormalization 

As anticipated in the general discussion of the formula ( |3-13j ) and seen explicitly in the 
previous subsection, the integral that formally embodies the three-point functions may be 
divergent. These divergences are understood and predicted by the results of holographic 
renormalization. Hence, it is possible to subtract the divergences and to obtain finite 
results. In this subsection, we shall illustrate this procedure using the integral for the 
three-point function (OSE). 

Before starting, we would like to emphasize also that the divergences are irrelevant for 
the physical three-particle scattering amplitudes. In fact, the integral (J4.34 ) involving the 



amputated bulk-to-boundary propagators are finite. This is easily understood by looking, 
e.g., at the bulk-to-boundary propagator K in the form ( 4.1 Q| ) and its amputated form 



(4.11). We have a factor (1 — u) in both formulae, which cancels a factor (1 — u) 1 in the 



integral measure. However, the hypergeometric function in ( |4. 10| ) goes like ln(l — u) for 
u — > 1 for generic a, whereas it is a polynomial on-shell. This means that, if the integral 
is divergent for a generic three-point function because of negative powers of (1 — u) in the 
integrand, in its amputated on-shell version we have a sufficient number of (1 — u) factors 
to cancel them. 

Let us start our illustration by summarizing the predictions of holographic renor- 



malization for the divergences. We shall use the Hamilton-Jacobi method [13], which is 
particularly simple to apply. Let us remind the reader that in the Hamilton-Jacobi method 
one first calculates the counterterms by solving a set of algebraic equations, and then one 



continues to obtain the exact one-point functions (3.5) 



Thus, in order to find the counterterms, we first expand the potential V{(p, a) about 



the fixed point (j) = a = 0. From (A. 3) we find 



a) = -3 - ~ (</> 2 + a 2 ) - 1 (0 4 - 3a 4 + 60V) + • • • . (4.41) 



After applying the procedure of solving the descent equations as explained in [13] we arrive 
at the following counterterms, 6 



S c . t .= j d 4 xVff|^ + i(0 2 + a 2 )-^-lne 



(4.42) 

__ (V^VV + V^VV) + - (4> 2 + a 2 ) R + - (r)R{ - ±R 2 ) J j . 

In order to find the exact one-point functions one writes down a series expansion for 
the fields. In the case of the field <fi we have 

4>{x,e) = 0(x)e 1/2 + 0(x)e 3/2 In e + (f){x)e 3/2 + ■■■ . (4.43) 

It is now clear that the logarithmic divergences of the three-point function integrals are 
understood as contributions to the logarithmic term (j) stemming from the interactions. 
Furthermore, it is possible to make an exact statement using holographic renormalization, 



6 The same counterterms were obtained in Q using the standard method of holographic renormalization. 
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because the term <f> is determined uniquely by the fact that the exact one-point function 
(O) is finite. For the counterterms ( 4.42| ) this yields [13] 



"3^ 



4 V ^"24^- 



(4.44) 



Focussing on (OSS), (f4.43j ) and ( [4.44 ) imply that the second term in ( [3.11 ) contains a 
term 

-^4>a 2 e 3/2 lne , 

from which in turn follows that the integral in ( |3.13| ) with the operator Xot,t, contains the 
logarithmic divergence 

8 - 8 

-(Alne = — lne . (4.45) 
3^ 

Similarly, it becomes also clear from the counterterms ( f4.42| ) that the integral for 
(OOO) is finite, the one for (TTT) is linearly divergent (the R counter term contains e*- 
to all orders), and the others that are not zero diverge logarithmically. These are precisely 
the divergences observed explicitly in subsection [I~2] . 

We shall now confirm the prediction ( 4.45| ) by calculating the divergence directly. It 
stems from the terms on the first line of ( 4.38| ). Substituting them into the integral ( 4.32| ) 
and considering only the leading behaviour of the bulk-to-boundary propagators ( [4.3j ) and 
Q4.15|) yields 



l-£ 



4u(2 - 3u) 8 u 2 pi - (p 2 + Pz) 



Pi 



1 — u 3 1 — u 
By virtue of momentum conservation this gives to leading order 



Vl - ud u VT 



(4.46) 



l-e 



du 



1 



VzJ i 



— = lne 



(4.47) 



which confirms (4.45). 

Having understood the origin of the divergence, one can proceed to remove it. This is 
done by subtracting the term from the integral that contributes to cj>, i.e., by adding 

8 1 



to the integrand. One could repeat the above discussion for all cases, but we shall be 
content with this illustration. 



5. Glueball Scattering Amplitudes and Discussion 

At this point we have at our disposal the necessary ingredients to calculate the scattering 
amplitudes for various three-particle processes of interest. In order to restore the proper 
physical dimensions, we should re-introduce the asymptotic AdS length scale L, defined 
by L 4 = 4ng s Na' 2 = Xta' 2 , where Xt = <7y M A^ is the 't Hooft coupling. This is simply 
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done by replacing p by pL everywhere, so that O(p) — > 0(pL) = 0{p)/L. {0{p) has length 
dimension —1 corresponding to 0(x) having dimension 3. The same relation holds for 
while Tj(p) has dimension 0.) Furthermore, the results for the correlation functions 
should be multiplied by a numerical factor [N 2 / (2ir 2 )] x (2-7r) 4 , where the factor N 2 /(2n 2 ) 
takes into account the 5-dimensional Newton constant, while the (27r) 4 stems from our 
convention for the (5-function in momentum space. We shall not do this explicitly, but 
only note that, after canonically normalizing the operators at the UV fixed point, which 
absorbs the factor TV 2 in the two-point functions, the three-point functions are suppressed 
by a factor 1/N. 

Before starting, however, we shall briefly summarize the particle spectrum found in 



subsection ^Jj in order to facilitate the ensuing discussion. 

The operator O, which is part of the M = 1 chiral anomaly multiplet, A = tr ( < 1 >2 ), 
creates the glueball states Ok of spin zero, which have masses 

(mi) 2 = 4k(k + 1) A; = 1,2, .... (5.1) 

The operator X is part of the "Lagrangian" multiplet S = tr (W 2 + • • • ) and creates the 
glueball states of spin zero, which have masses 

(mL) 2 = 4(k - l)(k + 2) A; = 1,2,.... (5.2) 

Finally, the operator Tj, which is part of the J\f = 1 supercurrent multiplet J a a = 
tr (W a Wa + •••), creates the spin-two glueballs 2\ with the masses 

(mL) 2 = 4(fc + l) 2 A: = 1,2,.... (5.3) 

The spin structure of the spin-two states can be described as follows. In spin space, 
the operator T % - is a symmetric, traceless and transversal 4x4 matrix, which forms a spin- 
two representation of the "little group" for massive particles, i.e., the subgroup of those 
Lorentz transformations that leave a given massive momentum four- vector p l invariant. In 
fact, it is instructive to choose as a basis a set of five matrices e r (p)p with r = 1,2, 3, 4, 5, 
satisfying the orthogonality relations 

£r(i?)i£ S (p)- = $rs , (5.4) 

in addition to the transversality and tracelessness conditions 

£r(p)}Pi = = e r {p)\ . (5.5) 
Using such a basis, the traceless transversal projector can be written in the form 

r=l 

Eqn. ( |5.6| ) can also be regarded as a completeness relation. The polarization matrices take 
a particularly simple form in the rest frame, where they are given by 

(5.7) 
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where the i r form a basis of symmetric and traceless 3x3 matrices. 

The spectrum of states includes also the superpartners of the states listed above, which 
have, obviously, the same masses. The spectrum of low-lying states is shown in Fig. [l|. 
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Figure 1: Spectra of states with mL < 15. 



We shall turn now to the scattering amplitudes, which, in an equivalent fashion, de- 
scribe the decay of a glueball into two other glueballs, or the creation of that glueball in a 
two-glueball collision. These processes are constrained by phase space such that the mass 
of the decaying (or the created) glueball must be at least as large as the sum of the masses 
of the decay products (or the colliding particles). It will turn out that the amplitudes for 
most processes, which are allowed by phase space, vanish. 

In order to simplify the discussion, we have chosen to consider the unpolarized ampli- 
tudes for processes involving the spin- two glueballs This is achieved by first projecting 
the generic amplitude containing indices i and j for each external spin-two state onto a 
particular polarization using the polarization matrices e r (p)*-, then taking the square and 
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finally summing over all polarizations r. For the amplitudes of TOO processes this implies 



a 

\M T oo\ 2 = \ e r{vifjM T oo\ e r (pi)fM T oo l k 

r=l 

= MTOoi^lflMTOOk , 

where we have used the completeness relation (fO|) . As Mtoo) is of the form 



(5.8) 



Mtoo) = (Hi %V2kv\ ) Mtoo 



we find from (|5 



\M T 00\' Z = [P2iP2kpip l 3^lTl) \MtOO 



ik 



{p\ +Pt+Pi~ 2p\p\ - 2p\p\ - 2plplf \Mtoo 



2Ap{ 



(5.9) 



The last line is easily established in the rest frame of the T particle. The unpolarized 
amplitude |7Wtes| 2 is obtained in the same fashion. 

Similarly, we might define the unpolarized amplitudes \Mott\ 2 and \Mttt\ 2 , which 
give rise to more cumbersome pre-factors involving the external momenta. We shall not 
provide their explicit expressions, because our numerical analysis will indicate that the 
amplitudes of all processes of these kinds, which are allowed by phase space, vanish. Hence, 
the overall factor is of no importance. 



The actual calculation of the scattering amplitudes ( 4.34 ) can be easily implemented 



on a computer. 7 The numerical results for the unpolarized glueball decay amplitudes for 
the glueball states with k < 10 are listed in appendix [B|. Here, we summarize our findings 
and the resulting decay channels. 

The states 0\ and Si turn out to be stable glueballs. For Si this is natural, because it 
is massless, but for 0\ it results from the fact that the only allowed process, 0\ — ► Ei + Ei, 
has a zero amplitude. In contrast, the process T\ — > Si + Si is allowed and occurs. 

A glueball Ok with k > 1 decays mainly into two S glueballs. It can decay into S^+Si, 
and into Sj + Sj such that i + j = k, but the latter processes are severely restricted by 
phase space. In fact, for k > 10, only the processes Ok — ► Sfc_ 2 + S 2 and Ok — ► S^_i + Si 
are allowed. Furthermore, Ok can decay into Oi + Tj such that i+j = k — 1, although these 
decay channels are much less probable. All other allowed processes have zero amplitudes. 
In particular, there are no decays of the form O — > T + T and O — > O + O. 

The decay of a glueball S^ with k > 1 must contain exactly one S glueball amongst 
the products. The main decay channels are into Sj + Tj such that i + j = k — 1, and into 
Ok-i + Si. The decay into Tk-i + Si also occurs, but has a much smaller probability. 
Again, all other allowed processes have vanishing amplitude. 

Finally, the glueballs Tk decay mainly into two S glueballs, the main channels being 
into Sj + Sj such that i + j = k, and into E& + Ei. Decays into Oi + Oj such that 



7 We have used MAXIMA. A script is available from the authors. 
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i + j = k are also possible, but less probable. As before, all other allowed decay processes 
have vanishing amplitude. 

Although our conclusions about the vanishing amplitudes of many decay channels stem 
from the numerical analysis of the decay amplitudes up to k = 20, we believe that there 
is a deeper reason, which is to be sought in some orthogonality relation of the Jacobi 
polynomials. We shall not try to give a more rigorous proof of these statements. 

A comparison of our results with similar data from lattice simulations of N = 1 SYM 
theory, when they become available, would be very interesting. 

The three-point function (OSS) has been calculated also by Bianchi and Marchetti 
Q. Although their formula differs from our ( |4.38| ) , it is possible to show, using integrations 
by parts and the equations of motion for the bulk-to-boundary propagators, that the two 
bulk integrals differ only by boundary terms. In the three-point function, these boundary 
terms would constitute contact terms, which we should drop, because none of us has done 
a reliable analysis of the contact terms. In the amplitudes, the boundary terms vanish 
because of factors of u and (1 — u). In fact, our numerical amplitudes agree completetly. 

In conclusion, the holographic analysis of the three-point functions for the operators 0, 
£ and Tj has yielded precise predictions for the glueball scattering amplitudes in the GPPZ 
flow. As the GPPZ flow shares some features with pure jV = 1 SYM theory (in particular 
confinement), it potentially sheds light on the IR dynamics of the latter. However, one 
should be cautious to draw too quick a conclusion. The GPPZ flow is a particular, unstable, 
case of a two-parameter family of Af = 1 holographic RG flow backgrounds describing the 
mass deformation of Af = 4 SYM theory H. In a generic member of this family both 
scalars considered in this paper, cj> and a, are active, and the latter describes a gaugino 
condensate, which is a necessary ingredient in the vacuum structure of Af = 1 SYM theory. 
However, all of these backgrounds are singluar, and it is unclear how to choose amongst 
the parameters the right values that describe an Af = 1 vacuum. By analogy with other 
gravity duals of Af = 1 SYM theory {e.g., the Maldacena-Nuhez solution |23||), one might 
argue that a truely 10-dimensional mechanism — unknown at present — will resolve the bulk 
singularities, thereby fixing the parameters and isolating the vacuum. 

Hence, one should regard the GPPZ flow as a toy model, whose qualitative features 
exist also in the theory with the true vacuum. These features include the existence of 
glueball states and the preferred glueball decay channels, although the numerics of the 
glueball masses and scattering amplitudes will be affected by the non-zero gluino condensate 
and the singularity resolution. It might, of course, have been better to consider a generic 
background of the two-parameter family of solutions in order to describe at least the effect 
of the gluino condensate. Unfortunately, there are two technical difficulties already at the 
linearized level making this problem much harder to tackle. First, the two active scalars 
present in these backgrounds couple to each other through the potential leading to a fourth 
order differential equation. Second, although the traceless transversal components of the 
metric decouple from all other fields, their equation of motion is not analytically solvable. 
Further progress on these issues is, therefore, very desirable. 



- 28 - 



Acknowledgments 



We would like to thank Massimo Bianchi for his collaboration in the early stage of this 
project and for helpful discussions, as well as Kostas Skenderis for sharing his experience 
in holographic renormalization with us. 

This work was supported in part by INFN, by MIUR (contract 2003-023852), by 
NATO (contract PST.CLG. 978785) and by the European Community's Human Poten- 
tial Programme (contracts HPRN-CT-2000-00122, HPRN-CT-2000-00148 and HPRN-CT- 
2000-00131). 

A. Useful Relations for the GPPZ Flow 

We summarize here a number of relations for the GPPZ background. For simplicity, we 
set the asymptotically AdS length scale to unity, i.e., L = 1. 

The potential V(4>, a) that gives rise to the GPPZ flow with (f> as an active scalar was 
found in It is given in terms of a superpotential W(<p, a) by 

v ( *,)-i(f)" + i(^V-^. 

where 



W{<P,a) = -- A 



2 V da J 

2(f)' 



cosh ( — — \ + cosh(2<j) 



Hence, we have 



coslr ( -^jQ +4cosh H^L ) rosh(2fr) - n>sh 2 ( 2a) + I 



(A.2) 
(A.3) 



The GPPZ background solution is 

a = 0, e 2A = e 2r - 1 . (A.4) 



l-e~ r 

From (A.4) we easily find the background source, 

$ = V3 . (A.5) 

For the GPPZ background (for a = 0), there are a number of identities that simplify 
the calculations with the potentials and its derivatives, namely 

4 Wl 4 /9 \ 

^=3^ + 1, "^ = 3^ + 2, V = -W[ 1 W-1), (A(j) 

V„ = 8W + 9 , V mt = SW^ , V c „ = . 

Finally, it is useful to introduce the variable 

u = 1 - e" 2r , (A.7) 
in terms of which the following relations hold, 



^ = 2(l-u), e- 2 ^ = i^, 
dr u 

2u u 



(Ai 
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B. Numerical Results for the Glueball Decay Amplitudes 



In this appendix we provide lists of all non-zero scattering amplitudes (|A4| 2 ) for the decays 
of the glueballs Ok, and with k < 10. Only the amplitudes for decay processes allowed 
by the phase space are given. The amplitudes involving T glueballs are the unpolarized 
ones as defined in Sec. ||. 
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Table 1: The unpolarized decay amplitudes for the glueballs with k < 10. 
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Table 2: The unpolarized decay amplitudes for the glueballs Ok with k < 10. 



- 31 - 







V 1 v T 1 V 1 

s 2 —> ±1 + Si 


U.U 


S 2 —> U\ + Si 


llo.z 


£3 


-> Ti + Si 


zy.ozo 


^3 - J 2 + ^1 


U.Uyl 


^3 —> ^2 + i-i 




s 4 


-> Ti + S 2 




V 1 T"' 1 V 1 

^4 — ► + ^1 


n /in 

0.49 


s 4 — »■ C 3 + Si 


TT 1 A O 

77.143 


s 4 


-> T 2 + Si 


A A flOO 










£5 


?i + S 3 


37.403 


V 1 ATI , V"1 

2^5 —> i 4 + 2ji 


1 r\ A TO 

1.0473 


S 5 -»■ O4 + Si 


/ion 

68.422 


£5 


— »■ T 2 + S 2 


r A 00 

oo.4oz 










£5 


T 3 + Si 


r OP. A 










s 6 


Ti + S 4 


OO A A f\ 

38.449 


^6 - ^ J 5 + ^1 


1.6275 


^6^05 + Si 


62.657 


s 6 


^ T 2 + S 3 


57.215 










s 6 


-> T 3 + S 2 


n A no 1 

64.081 










s 6 


-> T 4 + Si 


01 .016 










s 7 


Ti + S 5 


39.04 


S 7 -> T 6 + Si 


1 TOT 

2.1737 


S 7 -> 06 + Si 


r r n r 

58.595 


s 7 


-> T 2 + S 4 


59.152 










£7 


T 3 + S 3 


69.011 










£7 


-> T 4 + S 2 


70.982 










S 7 


-> T 5 + Ei 


p 1 i n 

61.349 










s 8 


Ti + S 6 


39.408 


\~i rri , \ ^ 

^ 8 —> + 7 + ^1 


2.6682 


S 8 -> C7 + Si 


r cr coo 

55.588 


s 8 


r 2 + s 5 


60.283 










s 8 


-► T 3 + S 4 


71.609 










s 8 


-► T 4 + S 3 


76.724 










s 8 


^ T 5 + E 2 


Tr* to ^ 

75.784 










s 8 


T 6 + Si 


c* A noo 

64.038 










Sg 


Ti + S 7 


39.652 


Sg — > T 8 + Si 


3.1094 


Sg -> Og + Si 


53.278 


£9 


-> T 2 + S 6 


61.003 










Sg 


^ T 3 + E 5 


73.161 










Sg 


-»■ T 4 + S 4 


OIO 

79.812 










Eg 


T 5 + S 3 


00 1 in 

82.119 










Sg 


^ T 6 + E 2 


1-7 c\ ono 

79.303 










Sg 


-► T 7 + Si 


nr* nop 

66.086 












^ T 4 + E 8 


39.822 


S10 — ► Tg + Si 


3.5011 




51.45 


^10 


— »■ T 2 + S 7 


61.49 










S10 


- T 3 + S 6 


/4.16( 










S10 


T 4 + S 5 


81.687 










S10 


T 5 + S 4 


85.577 










£10 


^ T 6 + E 3 


86.086 










£10 


-► T 7 + S 2 


81.987 










£10 


^ T 8 + Si 


67.698 











Table 3: The unpolarized decay amplitudes for the glueballs S& with k < 10. 
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